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Abstract 

The perturbative calculation of the lifetime of charged excitations in ultrarelativistic 
plasmas is plagued with infrared divergences which are not eliminated by the screen- 
ing corrections. The physical processes responsible for these divergences are the colli- 
sions involving the exchange of longwavelength, quasistatic, magnetic gluons (or pho- 
tons), which are not screened by plasma effects. In QED, the leading divergences can 
be resummed in a non-perturbative treatement based on a generalization of the Bloch- 
Nordsieck model at finite temperature. The resulting expression of the fermion propaga- 
tor is free of infrared problems, and exhibits a non- exponential damping at large times: 
Sfi(t) ~ exp{— aT t\nu p t}, where u> p = eT/3 is the plasma frequency and a = e 2 /An. 
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The perturbative calculation of the lifetime of charged excitations in ultrarelativistic 
plasmas is plagued with infrared divergences which are not eliminated by the screen- 
ing corrections. The physical processes responsible for these divergences are the colli- 
sions involving the exchange of longwavelength, quasistatic, magnetic gluons (or pho- 
tons), which are not screened by plasma effects. In QED, the leading divergences can 
be resummed in a non-perturbative treatement based on a generalization of the Bloch- 
Nordsieck model at finite temperature. The resulting expression of the fermion propaga- 
tor is free of infrared problems, and exhibits a non- exponential damping at large times: 
Sji(t) ~ exp{— aTt lnu; p t}, where u p = eT/3 is the plasma frequency and a = e 2 /47r. 

1. INTRODUCTION 

The study of the elementary excitations of ultrarelativistic plasmas, such as the quark- 
gluon plasma, has received much attention in the recent past. (See for recent reviews 
and more references.) The physical picture which emerges is that of a system with two 
types of degrees of freedom: i) the plasma quasiparticles, whose energy is of the order of 
the temperature T; ii) the collective excitations, whose typical energy is gT, where g is 
the gauge coupling, assumed to be small: g <C 1 (in QED, g = e is the electric charge). 
For this picture to make sense, however, it is important that the lifetime of the excitations 
be large compared to the typical period of the modes. 

Information about the lifetime is obtained from the retarded propagator. A usual 
expectation is that S#(i, p) decays exponentially in time, Sn(t, p) ~ e - j£; (p)* e -7(p)^ w here 
E(p) ~ T or gT is the average energy of the excitation, and ^(p) is the damping rate. 
Therefore, \Sji(t, p)| 2 ~ e _r ^* with T(p) = 2^j{p) ) which identifies the lifetime of the 
single particle excitation as r(p) = 1/T(p). The exponential decay may then be associated 
to a pole of the Fourier transform Sr(u,p), located at uj = E — 27. The quasiparticles 
are well defined if their lifetime r is much larger than the period ~ 1/E of the field 
oscillations, that is, if the damping rate 7 is small compared to the energy E. If this is 
the case, the respective damping rates can be computed from the imaginary part of the 
on-shell self-energy, £(u; = E(p),p). 

Previous calculations suggest that 7 ~ g 2 T for both the single-particle and the 
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collective excitations. In the weak coupling regime g <C 1, this is indeed small compared 
to the corresponding energies (of order T and gT, respectively), suggesting that the 
quasiparticles are well defined, and the collective modes are weakly damped. However, 
the computation of 7 in perturbation theory is plagued with infrared (IR) divergences, 
which casts doubt on the validity of these statements [3 — 9]. 

The first attempts to calculate the damping rates were made in the early 80 's. It was 
then found that, to one-loop order, the damping rate of the soft collective excitations in 
the hot QCD plasma was gauge-dependent, and could turn out negative in some gauges 
(see Ref. [[Kj for a survey of this problem). Decisive progress on this problem was made by 
Braaten and Pisarski || who identified the resummation needed to obtain the screening 
corrections in a gauge-invariant way (the resummation of the so called "hard thermal 
loops" (HTL)). Such screening corrections are sufficient to make IR- finite the transport 
cross-sections [TTyT2] ], and also the damping rates of excitations with zero momentum 
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At the same time, however, it has been remarked Jl that the HTL resummation is 
not sufficient to render finite the damping rates of excitations with non vanishing mo- 
menta. The remaining infrared divergences are due to collisions involving the exchange 
of longwavelength, quasistatic, magnetic photons (or gluons), which are not screened in 
the hard thermal loop approximation. Such divergences affect the computation of the 
damping rates of charged excitations (fermions and gluons), in both Abelian and non- 
Abelian gauge theories. Furthermore, the problem appears for both soft (p ~ gT) and 
hard (p ~ T) quasiparticles. In QCD this problem is generally avoided by the ad-hoc 
introduction of an IR cut-off ("magnetic screening mass") ~ g 2 T, which is expected to 
appear dynamically from gluon self-interactions 0. In QED, on the other hand, it is 
known that no magnetic screening can occur [14], so that the solution of the problem 
must lie somewhere else. 

In order to make the damping rate 7 finite, Lebedev and Smilga proposed a self- 
consistent computation of the damping rate, by including 7 also in internal propagators 
[|]]. However, the resulting self-energy is not analytic near the complex mass-shell, and 
the logarithmic divergence actually reappears when the discontinuity of the self-energy is 
evaluated at u> = E — ij |6|,[7j . More thorough investigations along the same lines led to 
the conclusion that the full propagator has actually no quasiparticle pole in the complex 
energy plane ||. These analyses left unanswered, however, the question of the large time 
behavior of the retarded propagator. 

As we have shown recently for the case of QED the answer to this question 

requires a non perturbative treatment, since infrared divergences occur in all orders 
of perturbation theory. We have identified the leading IR divergences in all orders, 
and solved exactly an effective theory which reproduces all these leading divergences. 
The resulting fermion propagator Sr(u) turns out to be analytic in the vicinity of the 
mass-shell. Moreover, for large times t ^> 1/gT, the Fourier transform Sn(t) does not 
show the usual exponential decay alluded to before, but the more complicated behavior 
<3R.(t) ~ e~ lEt exp{— aTt \nuj p t}, where a = g 2 /^ and u p ~ gT is the plasma frequency. 
This corresponds to a typical lifetime r _1 ~ g 2 T ln(l/ g), which is similar to the one 
provided by the perturbation theory with an IR cut-off of the order g 2 T. 
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Figure 1. The resummed one- loop self-energy 



2. THE INFRARED PROBLEM 

Let me briefly recall how the infrared problem occurs in the perturbative calculation of 
the damping rate 7. For simplicity, I consider an Abelian plasma, as described by QED, 
and compute the damping rate of a hard electron, with momentum p ~ T and energy 
E(p) = p. 

To leading order in g, and after the resummation of the screening corrections, 7 is 
obtained from the imaginary part of the effective one-loop self-energy in Fig. [l]. The 
blob on the photon line in this figure denotes the effective photon propagator in the HTL 
approximation, commonly denoted as *D^ u (q). In the Coulomb gauge, the only non-trivial 
components of *D flu (q) are the electric (or longitudinal) one *D 00 (q) = *Ai(q), and the 
magnetic (or transverse) one *Dij(q) = (5ij — qiqj)*A t (q), with 

where 11; and U t are the respective pieces of the photon polarisation tensor Physi- 
cally, the on-shell discontinuity of the diagram in Fig. [I] accounts for the scattering of the 
incoming electron (with four momentum p^ = (E(p),p)) off a thermal fermion (electron 
or positron), as mediated by a soft, dressed, virtual photon. (See Fig. |^.) 

The interaction rate corresponding to Figs. [I] or |2| is dominated by soft momentum 
transfers q <T. It is easily computed as [ 




where the upper cut-off q* distinguishes between soft and hard momenta: gT <C q* <T. 
Since the g-integral is dominated by IR momenta, its leading order value is actually 
independent of q*. 

The two terms within the parentheses in eq. (0) correspond to the exchange of an 
electric and of a magnetic photon respectively. For a bare (i.e., unscreened) photon, we 
have \Ai(q ,q)\ 2 = l/q A and \A t (qo,q)\ 2 = l/(<7o — ? 2 ) 2 > so ^ na ^ the g-integral in eq. @ 
shows a quadratic IR divergence: 

57T Jo q A 



Figure 2. Fermion-fermion elastic scattering in the Born approximation 



This divergence reflects the singular behaviour of the Rutherford cross-section for forward 
scattering. As well known, however, the quadratic divergence is removed by the screening 
corrections contained in the photon polarization tensor. We shall see below that the 
leading IR contribution comes from the domain qo <^ q ^ T, where we can use the 
approximate expressions |l]J^] (with u p = eT/3) 

IL(g < q ) ~ 3u 2 p = m 2 D , H t (q < q) ^ -i^-u 2 p ^. (4) 

We see that screening occurs in different ways in the electric and the magnetic sectors. In 
the electric sector, the familiar static Debye screening provides an IR cut-off ~ gT. 
Accordingly, the electric contribution to 7 is finite, and of the order ji ~ g A T 3 /m 2 D ~ g 2 T. 
Its exact value can be computed by numerical integration 0. In the magnetic sector, 
screening occurs only for nonzero frequency go Q- This comes from the imaginary part 



of the polarisation tensor, and can be associated to the Landau damping [16] of space-like 



photons (<2q < q 2 ). This "dynamical screening" is not sufficient to completely remove the 
IR divergence of 7$ , which is only reduced to a logarithmic one: 

9 4 T 3 [ q * A f q dgo 1 



24 Jo J- q 2tt g 4 + (37r^g /4g) 2 
9 *T nM g^^ 
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The unphysical lower cut-off \i has been introduced by hand, in order to regularize the 
IR divergence of the integral over q. The upper cut-off uj p ~ gT accounts approximately 
for the terms which have been neglected when going from the first to the second line of 
eq. (BJ). As long as we are interested only in the coefficient of the logarithm, the precise 
value of this cut-off is unimportant. The scale uo p however is uniquely determined by 
the physical process responsible for the existence of space like photons, i.e., the Landau 
damping. As we shall see later, this is the scale which fixes the long time behavior of the 
retarded propagator. 

The remaining IR divergence in eq. fl5|) is due to collisions involving the exchange of 
very soft (q — ► 0), quasistatic (q — ► 0) magnetic photons, which are not screened by 




Figure 3. A generic n-loop diagram (here, n = 6) for the self-energy in quenched QED. 



plasma effects. To see that, note that the IR contribution to 74 comes from momenta 
q <C gT, where \*A t (q , q)\ 2 is almost a delta function of q$: 



*MQo,q)\'' 



g 4 + (37ru;W4 g ) 2 
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(6) 



This is so because, as go ~~ > 0, the imaginary part of the polarisation tensor vanishes 
linearly (see the second equation (§)), a property which can be related to the behaviour 
of the phase space for the Landau damping processes. Since energy conservation requires 
go = qcos8, where 9 is the angle between the momentum of the virtual photon (q) and 
that of the incoming fermion (p), the magnetic photons which are responsible for the 
singularity are emitted, or absorbed, at nearly 90 degrees. 

3. A NON PERTURB ATIVE CALCULATION 

The IR divergence of the leading order calculation invites to a more thorough inves- 
tigation of the higher orders contributions to 7. Such an analysis reveals strong, 
power-like, infrared divergences, which signal the breakdown of the perturbation theory. 
(A similar breakdown occurs in the computation of the corrections to the non-Abelian 
Debye mass JT^.) To a given order in the loop expansion, the most singular contributions 
to 7 arise from self-energy diagrams of the kind illustrated in Fig. |3|. These diagrams 
have no internal fermion loops (quenched QED), and all the internal photons are of the 
magnetic type (the electric photons, being screened, give no IR divergences). Further- 
more, the leading divergences arise, in all orders, from the same kinematical regime as 
in the one loop calculation, namely from the regime where the internal photons are soft 
(g — > 0) and quasistatic (g — > 0). This is so because of the specific IR behaviour of the 
magnetic photon propagator, as illustrated in eq. (||). Physically, these divergences come 
from multiple magnetic collisions. 

This peculiar kinematical regime can be conveniently exploited in the imaginary time 
formalism (see, e.g., ||), where the internal photon lines carry only discrete (and purely 
imaginary) energies, of the form g = %to n = i2imT, with integer n (the so-called Mat- 
subara frequencies). The non-static modes with n^O are well separated from the static 
one go = by a gap of order T. I have argued before that the leading IR divergences 
come from the kinematical limit go — > 0. Correspondingly, it can be verified ]l5j that, 
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in the Matsubara formalism, all these divergences are concentrated in diagrams in which 
the photon lines are static, i.e., they carry zero Matsubara frequency. (To one loop order, 
this has been also notified in Refs. |J.) In what follows, we shall restrict ourselves to 
these diagrams, and try to compute their contribution to the fermion propagator near the 
mass-shell, in a non perturbative way. Note that, for these diagrams, all the loop integra- 
tions are three-dimensional (they run over the three- momenta of the internal photons), 
so that the associated IR divergences are those of a three-dimensional gauge theory. This 
clearly emphasizes the non perturbative character of the leading IR structure. 

As we shall see now, this "dimensional reduction" brings in simplifications which al- 



lows one to arrive at an explicit solution of the problem fll5| . The point is that three- 



dimensional quenched QED can be "exactly" solved in the Bloch-Nordsieck approximation 



Iql , which is the relevant approximation for the infrared structure of interest. Namely, 
since the incoming fermion is interacting only with very soft (q — > 0) static (go = 0) mag- 
netic photons, its trajectory is not significantly deviated by the successive collisions, and 
its spin state does not change. This is to say, we can ignore the spin degrees of freedom, 
which play no dynamical role, and we can assume the fermion to move along a straightline 
trajectory with constant velocity v (for the ultrarelativistic hard fermion, |v| = 1; more 
generally, for the soft excitations, v(p) = dE(p)/dp = v(p)p is the corresponding group 
velocity, with \v(p)\ < 1). Under these assumptions, the fermion propagator can be easily 
computed as JT3J 



S R (t,p) = i9(t)e- iE ^A(t), (7) 
where 

a / + n / 2rr, r v d ^ 1 1 -cost(v(p) ■ q) ) 

contains all the non-trivial time dependence. The integral in eq. @ is formally identical 
to that one would get in the Bloch-Nordsieck model in 3 dimensions. Note, however, the 
upper cut-off u p ~ gT, which occurs for the same reasons as in eq. (|5|). Namely, it reflects 
the dynamical cut-off at momenta ~ gT, as provided by the Landau damping. 

The integral over q has no infrared divergence, but one can verify that the expansion of 
A(t) in powers of g 2 generates the most singular pieces of the usual perturbative expansion 
for the self-energy ||15|| . Because our approximations preserve only the leading infrared 



behavior of the perturbation theory, eq. (^]) describes only the leading large-time behavior 
of A(t). Since the only energy scale in the momentum integral of eq. (|8|) is the upper 
cut-off, of order gT, the large-time regime is achieved for t ^> 1/gT. Note that, strictly 
speaking, eq. (|8|) holds only in the Feynman gauge. However, its leading large time 
behaviour — which is all what we can trust anyway ! - - is actually gauge independent 
and of the form (we set here a = g 2 and v (p) = 1 to simplify writing) 



A(w p t>l) ~ exp 1-aTt In u p t J. (9) 
A measure of the decay time r is given by 

- = «Tln^ p r = «Tfln^-lnln^%+ ...V (10) 

r p V aT aT J v ' 
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Since aT ~ gu p , r ~ l/(g 2 Tln(l/g)). This corresponds to a damping rate 7 ~ 1/r ~ 
g 2 Tln(l/g), similar to that obtained in a one loop calculation with an IR cut-off \i ~ g 2 T 
(cf. eq. (I)). 

However, contrary to what perturbation theory predicts, A(t) is decreasing faster than 
any exponential. It follows that the Fourier transform 

/oo roc 
dte- iut S R (t,p) = i dte it^-E(p)+i V ) A ( t ) ? (11) 
-00 Jo 

exists for any complex (and finite) u. Thus, the retarded propagator Sr(u) is an entire 
function, with sole singularity at Imw — > —00. The associated spectral density p(u,p) 
(proportional to the imaginary part of Sr(u, p)) retains the shape of a resonance strongly 
peaked around the perturbative mass-shell u = E(p), with a typical width of order 
~g 2 T\n(l/g) JT§. 



4. CONCLUSIONS 



The previous analysis, and, in particular, the last conclusion about the resonant shape 
of the spectral density, confirm that the quasiparticles are well defined, even if their mass- 
shell properties cannot be computed in perturbation theory. The infrared divergences 
occur because of the degeneracy between the mass shell of the charged particle and the 
threshold for the emission or the absorbtion of n (n > 1) static transverse photons. Note 
that the emitted photons are virtual, so, strictly speaking, the physical processes that we 
have in mind are the collisions between the charged excitation and the thermal particles, 
with the exchange of quasistatic magnetic photons. The resummation of these multiple 
collisions to all orders in g modifies the analytic structure of the fermion propagator and 
yields an unusual, non-exponential, damping in time. 

This result solves the IR problem of the damping rate in the case of QED. Since a similar 
problem occurs in QCD as well, it is natural to ask what is the relevance of the present 
solution for the non-Abelian plasma. It is generally argued — and also supported by 



lattice computations [19] — that the self-interactions of the chromomagnetic gluons may 
generate magnetic screening at the scale g 2 T (see and Refs. therein). As a crude model, 
we may include a screening mass /i ~ g 2 T in the magnetostatic propagator in the QED 
calculation. This amounts to replacing 1/g 2 — > l/(q 2 + /i 2 ) for the photon propagator in 
eq. (§). After this replacement, the latter equation provides, at very large times t ^ l/g 2 T, 
an exponential decay: A(t) ~ exp(— 7^), with 7 = aT\n(u p / fi) = aTha.(l/g). However, 
in the physically more interesting regime of intermediate times 1/gT <C t l/g 2 T, 
the behavior is governed uniquely by the plasma frequency, according to our result 
A(t) ~ exp(— aTt hxu p t). Thus, at least within this limited model, which is QED with 
a "magnetic mass", the time behavior in the physical regime remains controlled by the 
Bloch-Nordsieck mechanism. But, of course, this result gives no serious indication about 
the real situation in QCD, since it is unknown whether, in the present problem, the effects 
of the gluon self-interactions can be simply summarized in terms of a magnetic mass. 

To conclude, the results of Refs. |!5|Jl7l| suggest that the infrared divergences of the 
ultrarelativistic Abelian plasmas can be eliminated by soft photon resummations, a la 
Bloch-Nordsieck. For non-Abelian plasmas, on the other hand, much work remains to be 
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done, and this requires, in particular, the understanding of the non-perturbative sector of 
the magnetostatic gluons. 
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